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ABSTRACT 

Context. Multitransition CO observations of galaxy centers have revealed that significant fractions of the dense circumnuclear gas 

have high kinetic temperatures, which are hard to explain by pure photon excitation, but may be caused by dissipation of turbulent 

energy. 

Aims. We aim to determine to what extent mechanical heating should be taken into account while modelling PDRs. To this end, the 

effect of dissipated turbulence on the thermal and chemical properties of PDRs is explored. 

Methods. Clouds are modelled as ID semi-infinite slabs whose thermal and chemical equilibrium is solved for using the Leiden PDR- 

XDR code, where mechanical heating is added as a constant term throughout the cloud. An extensive parameter space in hydrogen 

gas density, FUV radiation field and mechanical heating rate is considered, covering almost all possible cases for the ISM relevant 

to the conditions that are encountered in galaxies. Effects of mechanical heating on the temperature profiles, column densities of CO 

and HtO and column density ratios of HNC, HCN and HCO^ are discussed. 

Results. In a steady-state treatment, mechanical heating seems to play an important role in determining the kinetic temperature 

of the gas in molecular clouds. Particularly in high-energy environments such as starburst galaxies and galaxy centers, model gas 

temperatures are underestimated by at least a factor of two if mechanical heating is ignored. The models also show that CO, HCN and 

H2O column densities increase as a function of mechanical heating. The HNC/HCN integrated column density ratio shows a decrease 

by a factor of at least two in high density regions with n ~ 10^ cm"', whereas that of HCN/HCO^ shows a strong dependence on 

mechanical heating for this same density range, with boosts of up to three orders of magnitude. 

Conclusions. The effects of mechanical heating cannot be ignored in studies of the molecular gas excitation whenever the ratio of the 

star formation rate to the gas density (SFR / «''^) is close to, or exceeds, 7 x 10"* M© yr"' cm*-' . If mechanical heating is not included, 

predicted column densities (such as those of CO) are underestimated, sometimes (as in the case of HCN and HCO^) even by a few 

orders of magnitude. As a lower bound to its importance, we determined that it has non-neghgible effects already when mechanical 

heating is as little as 1% of the UV heating in a PDR. 

Key words. GalaxiesiISM - (ISM:) photon-dominated region (PDR) - ISM:Turbulence - Physical data and processes:Mechanical 
Heating 
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en I 1 . Introduction Recently, it has become clear that PDR and XDR excitation 

^^ is not sufficient to fully e xplain observed mol ecular line emis- 

^^ sion ratios. For instance, iLoenen et al.l (l2008l) have found that 

T^_ , Radiation originating from the molecular gas in various galaxy the HCN, HNC, and HCO+ line ratios observed towards (ultra- 

^ ■ environments, including galaxy centers, provides information )iuminous infrared galaxies span a parameter space that can- 

• r-j ; about the physical state of these environments such as the gas ^^^ ^e reproduced by models in which FUV or X-ray radiation 

^ . mass and temperature, or the source of its excitation. The far- dominates. In their study they showed that an additional heat- 

H ■ ultraviolet radiation (FUV; 6.0 < E < 13.6 eV) emitted by jng mechanism is required, which they suggest to be mechan- 

. y: newly formed luminous stars and X-rays {E > 1 keV) pro- j^al heating caused by dissipating supernova shocks, injected 

duced by black hole accretion have very distinct effects on ^n large scales and cascading down throug h turbulent dissipa- 

the thermo-chemical balance of the gas. Strong FUV illumi- tion to the smallest scales. More recendy, lPaDadoDoulos"etan 

nation resuks in so-called nhoton-dominated regions (PDRs; ^i^ jj^s suggested that cosmic ray ionization rates in such 

iTielens & Hollenbach| (ll985|)), whi le X-ray irradia t ion cr eates galaxies may exceed those measured in the Milky Way by fac- 

X-ray dominated regions (XDRs; [Maloney et al.| (ll996|)). In ^rs of 10^ to 10^, hence significandy affect the ionization bal- 

PDRs and XDRs, the thermal and chemical structure is com- ance and thermal-chemical structure of interstellar clouds. The 

pletely determined by the radiation field, and direcdy reflects gff^ects o f high cosmic ray rates on ISM chemistry is studied in 

the energy balance of the interstellar gas. For example, in PDRs detail bv lMeiierink eFZI Muh (who co nsider mechanical h eat- 

the CO is generally much cooler (T~ 10-30 K) than in XDRs j^g as well only for selected models) and lBavet et al] (120111) . 

(r ~ 20 - 500 K), therefore XDRs generally exhibit CO emis- ___Z]Z~, 

sions to much higher J levels. Thus, observations of the CO lad- How ever, neither the models by iLoenen et al.l (l2008l) nor 

der allow us to discriminate between the different types of radi- those by iMeiierink et al.l ( 1201 ll) treat mechanical heating self- 

ation fields. consistently, as it is only implemented as an additional heating 
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term and not as feedback from hydrodynamics. The dynamical 
evolution of the interstellar medium in galaxies in general and 
close to a super-massive black hole in particular, has been stud- 



ied in 3D hydrodynamic simulations (IWada & Normani l2002t 
IWada & Tomisakal l2005h but these in turn lacked chemistry, 
and line emission maps of the various molecules were created 
by usin g constant abundances. A study bv lPerez-Beaupuits et al.l 
( 1201 lb i mproved on this by using che mical abundances obtained 
from the iMeiierink & SpaansI (l2005b XDR code. Nevertheless, 
these simulations still lack UV input and feedback from the 
hydrodynamics on the chemistry still needs to be done self- 
consistently. 

The dominant coolants in the ISM which are commonly ob- 
served are the [CI] 609 ^im, [CII] 158 jum and the [OI] 63 jum 
lines. These fine-structure lines are not sensitive to additional 
mechanical heati ng, when the clouds a re already illuminated 
by FUV radiation IMeiierink et al.l (|20II|) . This is in contrast to 
the molecular species that are formed in regions of the cloud, 
that are shielded from irradiation. Thus in this paper we con- 
sider molecular species as possible candidates for tracing re- 
gions which are dominated by mechanica l heating. Although a 
simpl e dust model is taken into account (iMeijerink & SpaansI 
l2005h . dust diagnostics are not discussed since these diagnostics 
are out of the scope of this paper We will explore the possi- 
ble tracers in a number of steps. In this paper, we will first es- 
tablish under which physical conditions mechanical heating is 
important. We will estimate what range of mechanical heating 
rates we may expect for the diff'erent densities in the ISM. We 
will then determine how this affects the thermo-chemical bal- 
ance in the models, and, by way of first application, discuss the 
effect on the integrated column densities of common molecular 
species such as CO, H2O, HNC, HCN and HCO+, and some of 
their ratios. In subsequent papers, we will calculate line emission 
maps that can be compared directly to the detailed observations 
obtained with (sub)millimeter telescopes, notably the Atacama 
Large Millimeter Array. Ultimately, our aim is to couple chem- 
istry and hydrodynamics in order to treat the relevant physical 
effects in a self-consistent manner. 



2. Methods 

We solve the equilibrium state of all model clouds by using 
an optimized versiorfl of the L eiden PDR-XDR code developed 
by IMeiierink & SpaansI (l2005l) (Details of the optimization and 
further improvements are mentioned in Sect. 12. 2) . In its present 
form, the code assumes a ID PDR geometry where various heat- 
ing and cooling processes are included in solving for the thermal 
and chemical balance. The radiative transfer for the cooling lines 
and the escape prob abilities are computed based on the work by 
Ide Jong et alJ (Il980h . We refer the reader to IMeiierink & SpaansI 
(i2005i) for further details on the thermal processes and chemical 
network used in solving for the steady state equilibria. 



2.1. Parameter space and reference models 

In order to reduce the number of free parameters, we consid- 
ered three metallicity levels, the benchmark being Z - Zq (so- 
lar metallicity) to which we added Z - Q.5Zq (characteristic of 
moderately metal-poor dwarf galaxies) and Z = 2Z0 (character- 
istic of spiral galaxy center conditions). Throughout the paper. 



our discussion concentrates on the Z - Zq modelo For the 
cosmic-ray ionization rate, we adopted the value 5 x 10^'^s"' 
used in [Meiierink & Spaans (2005) which is close to that in 
the solar neighborhood. The grids consisted of models span- 
ning a wide range of physical conditions, uniformly sampled 
(in logio scale) in ambient gas density 1 < n < lO^'cm"'', in- 
cident radiation field 0.32 < Gq < 10^ and mechanical heating 
lO"^'^ < Tmech < 10"'^ erg cm"^ s"'. In the absence of mechani- 
cal heating, « and Go dominate the characteristics of a PDR. The 
ambient gas density n is the total number density of hydrogen 
nuclei, which is set constant throughout the slab. The values con- 
sidered cover the full range of densities encountered in the dif- 
fuse interstellar medium to those of dense molecular cloud cores. 
The incident FUV rad i ation field Go is measured in so-called 
Habing units ( iHabind, 1 19691) corresponding to L6 x 10"^ erg 
cm"^ s '. The values considered here include radiation fields 
ranging from only a few times stronger than the relatively low 
solar neighborhood field to the intense fields irradiating clouds 
close to young OB associations; the upper limit corresponds to 
the radiation field intensity at a distance of about 0. 1 pc from an 
O star 

In constraining a physically relevant range of mechanical 
heating rates, we took a combination of values recovered from 
softened p a rticle hydrodynamics simulations (SPH hereafter) 
(iPelupessvL l2005l) . and added estimates of mechanical heating 
rates induced by turbulen ce caused by supern ova shocks that 
are dissipated in the ISM (iLoenen et all |2008|) . The SPH sim- 
ulations are of dwarf galaxies with n up to ~ 10^ cm"-* where 
the mechanical heating rates ranged from ~ lO"""^ to 10"^^ erg 
cm"^ s"' for such densities. When naively extrapolated to den- 
sities of 10^ cm""* 



r ates u i 

In the ILoenen et all ( 1200: 



to 10 ''^erg cm "* s 'are expected 
recipe, it is assumed that 77 = 10% 



'c 



of the supernova blast energy (£0 = 10^' erg per SN event) is 
absorbed throughout the continuous starburst region (a region 
which by definition sustains a considerably high SFR, which is 
steady, by consuming the available gas), whose size is taken to 
be Dfm - lOOpc . Eq.[T| which is a combination of Eq. 3 and 
4 in lLoenenet al.l (|2008|) . relates Fmech to the the supernova rate 
(SNR). 



SNR£o?7 

VsB VpDRnpDR 



(1) 



where Vsb is the volume of the starburst region, whereas Vpdr 
and «pDR are the volume and number density of PDRs in that 
starburst region respectively. They also assumed that the ab- 
sorbed energy is distributed evenly among the PDRs in this star- 
burst region. The number of the PDRs in the starburst region is 
estimated by assuming a gas density contrast of 10 (the inverse 
of the filling factor) between that of the PDRs and the ambient 
ISM of the starburst. Considering a Salpeter IMF for the stel- 
lar population, the star formation r ate (SFR) is related to the 
SNR via SNR/SFR = 0.0064 (see iDahlen & FranssonI (1 19991) 
for details). In computing the SNR, they as sumed that stars with 
M > 8 M0 end up as supernovae (see ILoenen et al. (l2008l) 
for mo re details on the estimates and assumptions). A iKroupal 
( l2001h IMF, over the mass range from O.IM© to 125M0 con- 
sidered, would result in a marginally higher SNR/SFR ~ O.OL 
This higher rate reflects a slightly enhanced mechanical heating 
rate. However, the difference in the resulting Fmech between both 
IMFs over the same mass range is ~ 40%, thus for simplicity we 
adopted the Salpeter IMF. 



The optimized code is called PDR-XDR- 1 . 1 



^ The analogous grids for the other two metallicities are presented in 
the Appendix 
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By way of example, a mechanical heating rate Fmech of 
2 X 10""'^ erg cm""' s"' is expected in a quiescent disc, such 
as that of the Milky Way with an SFR ~ 1 Mq yr"^ and 
an SNR of ~0.01 yr"' . In regions with active star formation 
(SFR~ 50 Mo yr ' and SNR~0.3 yr ') the mechanical heating 
rate is much higher, Fmech = 1x10"'^ erg cm""* s"' . For extreme 
starbursts with (SFR~ 1000 Mq yr"' and SNR~6.4 yr"') me- 
chanical heating rates F^ech = 2 x 10"'^ erg cm"-' s"' are possi- 
ble. Another source of mechanical energy input is the outflows of 
young stellar objects (YSOs). The outflow phase is short-lived, 
and as discussed by jLoenen et al.L l2008h a nearby dense cloud 
( « > 10^ cm"') is required for a significant amount of energy 
to be absorbed by the surrounding ISM (compared to the me- 
chanical heating due to supemovae). Hence energy input due to 
YSOs is ignored in this study. Based on a similar order of mag- 
nitude estimate, mechanical energy input due to stellar winds are 
also ignored, since t hey would contribute up to 6% to the total 
mechanical heating dPelupessvi 120051 : iLeitherer et al.L Il999h . 

It is of interest to determine the effects of mechanical heat- 
ing by comparing the results obtained from models with and 
without a mechanical heating term. To this end, we have se- 
lected from the grids calculated for analysis and discussion in 
this pa per the same reference mod els that were studied in de- 
tail by iMeijerink & SpaansI (l2005l) without taking into account 
mechanical heating. These models are summarized in Table[T] 



Model Name 


log 


n (cm ') 


log Go 


Ml 




3.0 


3.0 


M2 




3.0 


5.0 


M3 




5.5 


3.0 


M4 




5.5 


5.0 



Table 1. Para meters of the reference models. These are identical 
to the ones in iMeiierink & SpaansI ( |2005|) . Ml and M2 corre- 
spond to low density clouds in galactic centers in the presence 
of a starburst event, whereas M3 and M4 correspond to much 
denser clouds where excitation of high density gas tracers such 
as HCN is possible. 



In this first exploration, Fmech is added per unit volume. In 
doing so for a certain grid, each model in the n - Go parameter 
space has an extra F^ech added to its heating budget. Another 
choice could have been the addition of Fmech per unit mass, 
where for a certain grid the amount of added mechanical heating 
would be proportional to the gas density of the model in the grid. 
For simplicity we considered adding mechanical heating per unit 
volume. 



2.2. Numerical schemes 

Finding the thermal and chemical equilibrium for each slab is an 
(A^+ l)-dimensional numerical root finding problem. A^ being the 
number of species (atomic and molecular species and electrons) 
in the chemical network in addition to the thermal balance. It 
involves solving the set of non-linear equations Eqs.|2]and|3] 



dn(T)/dt = 

r(r, n) - A(r, n, o = o 



(2) 
(3) 



by iTielens & HagenI (Il982h in using the multi-dimensional 
Newton-Raphson method. The thermal balance is solved for in 
the second equation, where F and A are the total heating and 
cooling rates. These two terms are a ID function of n (used 
as input from the previous equation) and T. This process is re- 
peated iteratively until convergence with a tolerance of 10"^ for 
the chemical balance and 10"' for the thermal balance. Two ma- 
jor improvements lead to a speed-up by a factor of ~ 30 over 
the original PDR-XDR code. These are (a) tweaking the LU 
decomposition (optimizing compiler flags and matrix elements 
storage scheme) used in the multi-dimensional Newton-Raphson 
root finding, and (b) using in the thermal balance root finding the 
secant method (with a convergence order of ~ 1 .6) instead of the 
bisection method which has a linear convergence rate. 

The distance from the surface of the cloud is measured in 
terms of the visual extinction (Ay) due to interstellar dust. It is 
related to the tot al column density o f hydrogen nuclei N{li) ac- 
cording to Eq.l4l( lBohlin et al.Lll978b . 



Av^5.34x 10"^^A^HZmagcm"^ 



(4) 



Once the equilibrium state of the surface slab (with zero thick- 
ness) is solved for, the width of the remaining slabs is chosen 
adaptively while resolving the transition zone which has large 
temperature gradients compared to the small temperature gradi- 
ents at the surface of the slab and at high depths. This is also 
one of the changes over the original implementation. In choos- 
ing the slab width adaptively, the maximum allowed relative dif- 
ference in the temperature between consecutive slabs was set to 
5%. Smaller values were tested as well, where no significant dif- 
ference was noticed. 

It is assumed that the mechanical heating due to turbulence 
is absorbed by the ISM at all scales. Hence for simplicity Fn,ech 
is added uniformly throughout the cloud. 

Molecular clouds may ha ve a visual extinction up to 200 
mag (iTielens & HagenL Il982h . Beyond Ay ~ 5 mag most of 
the species become molecular whereas for Ay > 10 mag all 
abundances are almost constant since physical conditions do not 
change anymore. The maximum depth, we allowed for, in all 
the models, was Ay - 20 mag. This corresponds to A^(H) = 
7.5 X 10^2, 3.7 X 10^2 and 1.8 x lO^^^m-z for metaflicities of 
Z =0.5, I and 2Zq respectively. This enables us to compare 
clouds of different metallicities up to the same column density. 
In doing so, a fixed A^(H) = 1.8 x 10^^ cm"- (the maximum 
A^(H) for Z =2Zo) was adopted in the illustrations where column 
densities of molecular species were used. 

For each value of mechanical heating and metallicity, 25 x25 
models equally spaced in logio« and logjoGo were calculated, 
up to a total of ~16000 models. When there is no or very little 
mechanical heating, the temperature may drop below 10 K at a 
certain depth in the cloud. Since the code is designed to work for 
10 < r < 10"* K|j, for slab elements where the lower limit is 
reached we solved the chemical network for a fixed temperature 
of 10 K. In making the convergence more robust, we built a small 



The first equation (chemical balance) is the rate of change 
of the density of the species at a certain temperature T. We 
solved for this vector equation of length A^ following the recipe 



^ The lower bound in the temperature is set to 10 K since man y re- 
action constants become inaccurate below that value JRQlligLllOllh . He 
and Fe+ cooling become important above 10"* K which make the thermal 
balance inaccurate since these are not included in the current version 
(in addition to other dielectronic recombination processes). But below 
10"* K, the contribution of these processes is negligible to the total cool- 
ing. Also for T > W* K, the heating efficiency expression becomes 
inaccurate. However, not many models exceed 10"* K, but we allowed 
that since it is needed to make the models converge as the solution is 
progressed into the cloud. 
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database of guess values (used as input in the root finding) for the 
surface slab. This enabled us to explore regions in the parameters 
where the original implementation failed. 

In running the models, an improved version of the PDR- 
XDR code was incorporated into the Astrophysical Multi- 
purpose Software Environment (AMUSE) framework. The 
AMUSE package allows astrophysical codes from different do- 
mains to be combined to co n duct n umerical experiments (see 
iPelupessv & Portegies Zwarg (1201 2h for a more complete de- 
scription). It is a development of the M ulti-physics and Multi- 
scale Software Environment (MUSE, iPortegies Zwart et al.l 
|2009|) and is freely available for downloaqj. (The PDR code it- 
self will be made available in the future). The interface to the 
Meijerink PDR-XDR code takes the main parameters Go, n, 
Tmech as iuput and calculates the PDR equilibrium properties. 
This will allow the code to eventually be used as an equilibrium 
sub-grid model for e.g. galaxy scaled hydrodynamic simulations. 
In the present work this feature is not used, only the capability 
to run many models in parallel is. 

3. Results 

In this section we present the comparisons for thermal and chem- 
ical diagnostics of the grids modeled with and without mechani- 
cal heating. We start by discussing the temperatures at the cloud 
surface, and at the maximum cloud depth considered since for 
a given model, the gas temperature to a large extent determines 
the chemical composition. Following this, we will discuss the 
statistical properties of the CO, H2O, HCN, HNC and HCO+ 
molecules which are important diagnostics for moderate and 
dense gas. Throughout this section, we use the reference models 
(Ml to M4) to illustrate the impact of mechanical heating on the 
properties of the PDRs as a function of depth. 



te mperature dependen ce on metallicity has also been observed 
bv lRollig et al.l (l2006h . It is worth noting that in the lower right- 
hand corner of the grid (n > 10^^ cm"^. Go < 5) three equi- 
libria were found, one unstable and two stable ones. The stable 
equilibria correspond to a low (~ 50 K) and a high temperature 
(~ 300 K) state with the unstable one in-between. Here we have 
shown the low-temperature equilibrium. In any case, we note 
that this part of the parameter space is not of physical relevance 
for systems that interest us here.. 

The dominant heating term at the cloud surface is photo- 
electric heating, where electrons ejected from the dust grains 
collide elastically with the gas and heat it. The typical effi- 
ciency (e) of thi s process is quite low 0.5 to 1% with a max- 
imum of ~ 5% (iBlack & Dalgarng Il977h . The photo-electric 
heating rate (F photo) at the surface is given by the expression 
(iBakes & Tiele ns. 1994): 



photo 



lO^^'^eGon erg cm^^'s"' 



(5) 



where e depends on electron density {lie), G p and gas temper- 
ature . The expression for e is given in Eq-|6] dBakes & TielensL 
11994 . 



3.65 X 10-2 (r/lO"*)'" 
1 + 4 X 10-3 (Gori/2/«/" ^ 1 + 2 X 10-4 (Gori/2/«,) 



4.87 X 10-2 



(6) 



Hence the logn) of Eq. |5]is linear in logjo n and logjo Go when- 
ever e is constant. Fig. |2] shows the total heating in the n - Go 
plane. In the top part of the plot, the contours are almost straight 
lines. This is not surprising, Fphoto accounts in most cases for 
more than 95% of the total heating budget. To illustrate this, we 
look at the logarithm of Eq.|5] which in expanded form becomes 



3.1. Gas temperature 

In the absence of other heating terms (such as due to mechanical 
heating) PDRs are characterized by rapid decreases in temper- 
ature, from ~1000 K at the surface down to ~ 10 K in regions 
where Ay exceeds 10 mag. Since mechanical heating may be 
expected to considerably change the temperature structure, we 
will first explore this. 



3.1 .1 . Temperature at the surface 

At the surface we only need to solve for the equilibrium of a sin- 
gle slab (with zero thickness). Because this can be done much 
faster than in the case of a full slab, a higher resolution is pos- 
sible for the surface grids. These surface grids therefore con- 
sisted of 2500 (50 X 50) models uniformly sampled in logi()« 
and logjQ Go. In Fig.[T]the surface temperature, without mechan- 
ical heating (Fn,ech = Oerg cm-^ s-'), is shown for three differ- 
ent metallicities. The temperature c ontours are in good ag ree- 
ment with the profiles recovered by iKaufman et al.l (1 19991) and 
[Meijerink et al. (2007)(which uses the original implementation). 
We note that metallicity has little effect on the overall topology 
of the contour plots, except for regions where Go > 10"^, espe- 
cially where n > 10"* cm-^. The temperature almost doubles as 
Z is increased from Z© to 2Z0, with no noticeable change be- 
tween Z = O.5Z0 and Z - Zq. This is because cooling scales 
as Z whereas the dominant heating term scales almost as Z for 
low n but increases to ~ Z^ as n increases to 10^ cm--'. This 




log, Jn/cm" ] 

Fig. 2. Total heating at the cloud surface without mechanical 
heating. The main characteristic of this plot is that above the 
diagonal, the contours are parallel straight lines; whereas below 
the diagonal (Go = n line), there is a break in the slope, indicat- 
ing that photo-electric heating is no longer the dominant mech- 
anism in that region. The minimum total heating is of the order 
of 10-2" 
cm--' s 



erg cm 
-3 "-' at the top right corner 



s ' at the lower left corner and 10 '^2 erg 



www . amusecode . org 



logio Tphoto = -24 + logio e + logjo Go + logi,, n 



(7) 
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Fig. 1. Temperature profile (without mechanical heating) at the cloud surface for different metallicities. Generally, the temperature 
rises as a function of Go and constant n, since the dominant heating term is the photo-electric heating, which is prop ortional to Gq. 
H owever the behavior fo r constant Go as a function of n is not as trivial. See Section 3.1 of iKaufman et al.l d 19991) and Section 3 
of iMeiier ink et al.' (2007) for more discussion of the topology of the contours. This version of the code covers a larger parameter 
range than Meiierink et al. (2007) . Solid lines with filled squares highlight the temperature contours for 100, 500 and 1500 K . The 
most significant change when comparing the three panels is for the 1500 K contour, where for Z = 2Zq higher temperatures are 
attained ( compared to the lower metallicities) in the same part of the parameter space ( « > 10^ cm"-^ , Go > 10^ ) . 



which is linear in both logjgGo and logj^n. This means that 
the heating efficiency is almost constant wherever the contours 
are straight lines. This is not true in the regions below the 
break in the contour lines which occurs in the region under the 
Gq - n line in Fig.|2] In that part of parameter space, photo- 
electric heating is complemented by H2 photo-dissociation heat- 
ing. However, this is a small part of the total parameter space, so 
that we are justified in using photo-electric heating as a good ap- 
proximation for the total heating. Hence, we expect total heating 
to have a linear dependence on logjQ n for constant values of Go. 
This is achieved by fitting horizontal cuts of Ftotai in Fig.|2]as a 
function of « for different values of Go. 

From the fits, we could see that the minimum surface heating 
for the lowest density gas considered (« = 1 cm"^) ranges from 
jQ-25.6gj.g^jj^-3 g-i (for low Go) to lO^^^'^ergcm"-' s"' (for high 
Go). In contrast, the maximum surface heating for n - 10^ cm"-', 
is 10"'^'^ erg cm""* s"'. We also observe that the surface heating 
scales as «'-' for Go < 10- and n'-^^ for Go > lOl 

The equilibrium temperature is expected to vary depending 
on the amount of extra heating introduced via Fmech- Following 
this simple assumption, we expect to have significant changes 
in the equilibrium temperature when Fmech is comparable to the 
total surface heating, with a varying impact on different parts of 
the parameter space. In order to determine the zone in the pa- 
rameter space where Fmech alters the equilibrium state, we solve 
for n in Eq.|2]by equating it to Fmech- In other words, we solve for 
n in Ftotai,surface(n) = Fmech for the range in Fmech we have consid- 
ered in our parameter space. These solutions, n^., listed in Table|2] 
mark the lines in the Go - « plane where mechanical heating is 
the same as the total surface heating. To the left of these lines, 

A mech '^ A surface- 

In Fig.[3]the equilibrium temperature of the surface slab is shown 
for the range of Fmech that we consider Regions, that are first 
affected by the mechanical heating input, are the ones where 
Fmech is Comparable to the surface heating without mechanical 
input (see Fig-|2]i. This behavior can be traced by following the 
location of the dashed vertical lines (where Fmech = Ftotai.surface) 
and comparing the values of Fmech to the total surface heating 
contours in Fig.|3] It is clear that whenever the value of Fmech is 
comparable to the total surface heating, the temperature contours 
become almost vertical in the n - Go plane in Fig. [3] The loca- 
tions where Fmech = Ftotai.surface are marked by the dashed vertical 
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Table 2. Estimated values of gas density, where clouds with 
lower densities (listed in the bottom row), the shape and equi- 
librium temperature of the contours of Fig. [T] are significantly 
effected by Fmech (top row). For Fmech < 10"'^ erg cm""* s"', the 
fitting parameters for logjg Go = 2 are used. For higher mechan- 
ical heating rates, those for logjyGo = 5 are used since in that 
region of the parameter space Fmech is more relevant to high Go. 



lines where n = ric ( see Table|2]i. For instance, in the first panel 
in Fig. |3] (Fmech - 10"^"^ erg cm"-' s"'), regions where temper- 
atures increase significantly are when Ftotai.surface < lO"^'' erg 
cm"^ s"'(see Fig. [3] and compare first panel in Fig.|3]to the mid- 
dle panel in Fig-[T]i. Also we notice that the temperature increases 
from ~ 3000 K (for Go ~ 10^) to ~ 10000 K at the density 
« = lie, whereas it increases from ~ 5000 K (for Go ~ 10^-") to 
~ 7000 K at the same density. This can be explained by the fact 
that heating at the surface is mainly due to photo-electric heat- 
ing, thus as Go increases (for a fixed n), the amount total heating 
increases. Hence a fixed amount of mechanical heating would 
contribute less to the total heating budget, leading to a lower in- 
crease in the surface temperature. Subsequent panels in Fig.|3] 
which correspond to higher amounts of Fmech can be interpreted 
in the same fashion. The main feature of these panels is the shift 
in the position of the line marking n - ric towards higher den- 
sities. In the last panel where Fmech = 10"'^ erg cm""* s"' , the 
whole parameter space is dominated by mechanical heating. 

Regions with a SFR of ~ 1 M© yr 'are dominated by me- 
chanical heating up to densities of ~ lO"* cm""* . In Extreme star- 
burst, where SFR of ~ 1000 Mq yr^'are possible, all clouds with 



« < 10 cm become dominated by mechanical heating. In 
both cases, surface temperatures of > 5000 K are reached, and 
the temperature contours become independent of FUV heating. 



3.1 .2. Temperatures deep into the cloud 

Now we turn our attention to the molecular region. The solution 
was terminated well beyond Ay ~ I mag (where mainly atomic 
and radicals exist) at Av = 20 mag (where everything is molec- 
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Fig. 3. Surface temperature (for Z = Z©) for different values of mechanical heating. Each panel shows the surface temperature 
contours for a diflFerent value of Tmech mentioned at the top of each panel. In probing the effect of mechanical heating, each panel 
should be compared to the middle panel of Fig. [1] which is also for Z = Zq. Different parts of the parameter space in n and Go are 



affected differently, depending mainly on the gas density. For Fmech < 10"^' erg cm"^^ s"' , regions in density with n < lO"* cm^^'have 
their temperatures increase significantly. The vertical dashed red lines mark densities n - «(. where F^ech - rtotai.surface- The location 
of these lines moves from very low densities {n ~ 10*^.5 cm""*) to very high densities ( « > 10^ cm""*) as Fmech increases from 10"""* 
to 10"'^ erg cm"-' s"', sweeping the parameter space in density (for a certain grid) from left to right. In general the contour lines 
at « = He become vertical, indicating that the surface temperature becomes independent of Gq. In subsequent plot, those lines are 
plotted only for Z -Zq grids, since the fits from which n — tic were derived were done only for solar metallicity. 



ular). As the FUV radiation penetrates the cloud, it is attenuated 
by dust absorption. Thus, the photo-electric heating that was the 
dominant at the surface, rapidly becomes less important and al- 
most vanishes because (a) Fphoto has a exp {—Ay) dependence on 
Ay , and (b) the heating efficiency e becomes very small as elec- 
tron abundances decrease significantly by recombination. Thus, 
in the absence of Fmech, heating by cosmic rays and dust (espe- 
cially at high Go) become the dominant ones in the molecular 
region. 

In Fig.|4]we show contours of the equilibrium temperature 
at Ay - 20. In the absence of any mechanical heating, and in 
contrast to the situation for the surface slab, the topology of the 
contours is quite simple. The highest temperatures, with a maxi- 
mum of ~ 80 K, are in the upper right corner at n - 10^ cm"^ and 
Gq = 10^. Temperatures decrease monotonously along the di- 
agonal until they reach 10 K, which is the minimum tempera- 



ture we allowed for The main characteristic of a PDR is that 
the surface is heated to high temperatures, which is a thin layer 
Beyond this layer, in the transition zone, the temperature drops 
sharply. Although the temperature profile in Fig.^Jis simple, it is 
counter-intuitive. We might expect the temperature to decrease 
as the density increases at a fixed Go, but the opposite occurs. 
We explain this by the fact that as n increases, the coupling be- 
tween gas and dust also increases. Dust temperature depends on 
Go and decreases as a function of Ay- Nevertheless, for very high 
Go > 10^^, the decreasing dust temperatures are still high enough 
to heat the gas (~ 30 K). 

The temperature profile throughout the cloud is shown in 
Fig-Hfor the selected models Ml to M4. The temperature gen- 
erally decreases as a function of Ay. For some models which 
have high mechanical heating rates, the temperature increases 
and then saturates at higher values than those obtained at for 
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Fig. 4. Temperature al Ay - 20 without mechanical heating. 
Higher equilibrium temperatures are attained for dense molecu- 
lar clouds (top right comer) with decreasing equilibrium temper- 
ature towards low n and Go- Below the 20 K contour, all tempera- 
tures are set to 10 K - the minimum allowed in the code. The red 
dashed lines mark the densities n - n^ where Fmech = ^ total, surf ace 
from Table|2] 



the surface. As for the models whose temperature decreases, we 
see that as the density decreases, more specifically as the ratio 
Go/« increases, the point of saturation of the temperature gets 
shifted further away from the surface of the cloud [Ay -0 mag). 
For instance, when looking at clouds illuminated with a high 
FUV flux of Go = 10^, the saturation point gets shifted from 
Ay ~ 5 mag for n = lO-' cm"-' (low density gas) to Av ~ 2 mag 
for « = 10^^ cm""* (high density gas, compare left panels of 
Fig.|5j. A similar behavior is seen for the right panels which 
have a higher FUV flux with Go - 10^. However, we note that 
the saturation point remains almost constant for a fixed Go/«, 
which can be seen in comparing Ml and M4, where Go/« are of 
the same order of magnitude. The temperature saturation point 
is also related to the H/H2 and C^/C/CO transition zones. For 
example, along the lines (in n - Go parameter space) where 
Go/« - 0.003, 1 and 300, the location of the transition zones 
in the chemical species are almost unchanged. 

In a PDR, the total heating at the surface it is dominated 



and M2, and by at least one order of magnitude for the high 
density models M3 and M4. In fact, it suffices to have a me- 
chanical heating rate 10 to 100 times smaller than the total heat- 
ing at the cloud surface for gas temperatures to increase by a 
factor of ~2 (for high densities) and up to a factor of ~10 (for 
low densities). Thus even relatively small mechanical heating 
rates will substantially modify the physical parameters deduced 
from a pure PDR model. In a continuous starburst with an SFR 
- 1 M0 yr"' (which corresponds to Fmech ~ 2 x 10"^° erg cm"^ 
s ' using Eq. [TJ, the mechanical luminosity (heating) is about 
1% of the bolometric luminosity (iLeitherer et al.l [l999l) . Which 
implies that the ratio of mechanical luminosity to UV luminosity 
(where the emission of star forming regions peak) is larger than 
1 % . Especially if we consider higher a SFR ~ IOMq yr"' which 
translates to a Fmech ~ 2 x 10"'^ erg cm"-' s"' . Assuming that 
half the UV radiation is absorbed by the ISM of the galaxy, we 
would still be in the range where Fmech is at least 1 % of the sur- 
face heating (which is due to UV radiation), thus significantly 
altering the equilibrium temperatures of PDRs. 
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Fig. 5. Gas temperature as a function of the visual extinction for 
Ftotai, Av=o = rsiirface, and decreases as a func- selected models. In the top two panels, for F^ech > lO^^^erg 

cm"^ s"' (green and orange curves) the transition from H to H2 
does not occur For the high density models the transition fails 
only for M4 when Fmech > 10"'^ erg cm"^ s"'. 



by Fphoto, Av=o 

tion of increasing Ay- Also total heating generally decreases as 
a function of Ay. In Section. [3.1. H it was demonstrated that me- 
chanical heating increases the surface temperatures significantly 
whenever F^ech - Fsmface ■ In such situations, since mechanical 
heating is added uniformly throughout the PDR, we expect it to 
have a stronger impact (on the physical and chemical properties) 
at high depths compared to those at the surface. Hence Fmech 
should be the most dominant heating term beyond the surface, 
even more dominant than heating due to cosmic rays which at- 
tains a maximum of 10"^'-^ erg cm"-' s"' for « = 10^ cm"'. In 
Fig. lA. II we show the corresponding parameter grids as a func- 
tion of Fmech; the thermal profile in these grids has almost no 
dependence on Go and depends solely on n. 

In order to illustrate the degree by which Fmech increases tem- 
peratures throughout the cloud, we look at the curves in Fig. [5] 
corresponding to Fmech ^ 0. We note that for densities where 
Fmech = Fsurface, Hsted in Table|2] gas temperatures increase by 
almost two orders of magnitude for the low density models Ml 



3.2. Molecular gas tracers 

Now that we have determined the degree to which mechanical 
heating alters the thermal properties of PDRs, we examine its ef- 
fect on the chemical properties. Molecular clouds consist mainly 
of molecular hydrogen which is not directly observable except 
under very special circumstances. Thus, we must depend on the 
(sub)millimeter line emission from various transitions of much 
less abundant tracer molecules. Each molecular line represents 
its own critical density and temperature. 

In the following sections, we will consider the effect of me- 
chanical heating on a few species commonly used as tracers. 
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These are CO for relatively low gas densities and HCO^, HCN, 
and HNC for relatively high densities, in addition to H2O which 
has been also recently observed in very dense gas environments 
« ~ 10'' cm^^ 

3.2.1. CO 

Properties of molecular H2 gas are generally studied through ob- 
servations of CO, which is the next most abundant species. The 
critical density (nco) for the ground-state 7=1-0 '"CO transi- 
tion is yicr ~ 2000 cm"^. This is in the middle section of range 
of densities n considered in our grids. Thus, the response of CO 
to mechanical heating may have important consequences for es- 
timates of the amount of molecular gas and, for instance, may 
also explain the high gas temperature of 100-150 K which have 
been estimated for molecula r gas in the c enters of galaxies with 
densities of 100-1000 cm^^ (llsraelLl2009l) . 

In Fig.|6]we show the abundance of CO as a function of depth 
inside the cloud. For the models with densities close to n^, (Ml 
and M2), CO abundances decrease sharply for Tmech > 10"^ erg 
cm"-' s"' and almost disappears when Fmech > 10"^'^ erg cm"-' 
s"', where jcco < 10"'^ caused by very high temperatures 
T > 10000 K (abundances of species relative to the total number 
of hydrogen nuclei are denoted by xi = ni/{nu + 2«h,)). At such 
high temperatures, no molecules, including CO, are formed. This 
is due to the very low abundance of H2 which drops from ~ 0.5 
to ~ 10"^ beyond the radical region. This drop breaks the chem- 
ical network which leads to the formation of CO. In addition, 
the higher abundances of the ionic species H3 and He+ enhance 
the destruction of CO through the reactions H3 + CO —> HCO^ 
+ H2 and He+ H- CO ^ O + C+ + He. On the other hand, for 
Tmech ^ 10"'^ erg cm""* s"' the first obvious effect of mechan- 
ical heating on CO is the shift of the C/CO transition zone to- 
wards the surface of the slab. This shift is caused by the enhanced 
abundance of O2 and OH and higher temperatures in the radical 
region (see Fig.|6]l. 

In the high density models M3 and M4, the abundance of CO 
decreases by a factor of two as r^ech increases from 10""^ erg 
cm"-' s"' to 10"'^ erg cm"-' s"'. As F^ech is pushed further to 
10"'^ erg cm"-' s"' xco decreases from 10"'' to 10"'^ for M3 
and becomes under-abundant for M4 {xco < 10"'^). This be- 
haviour can be explained by the enhanced abundance of the 
ionised species He^ and H3, which destroy CO through ion- 
neutral reactions. Also, the significant drop in abundance of H 
and HCO (by three orders of magnitude) reduce the formation 
of CO through the reactions H + HCO ^ CO + H2 and O + C2 
^ CO + C. 

In Figs.r7land lA.2l the column density of CO is shown with- 
out and with increasing values of mechanical heating respec- 
tively. The general trend is an increase in N{CO) whenever the 
added F^ech is comparable to the surface heating of a PDR with 
no mechanical heating. This occurs at densities (n^) marked by 
vertical dashed lines in Fig. lA.2l For higher densities in the grids 
(to the right of these lines) F^ech < Tsurface and the opposite is 
true for lower densities (to the left the lines). At densities a fac- 
tor of 3 higher than /if (a distance of 0.5, in logm scale, to the 
right of «f ) where logj,, « = logjo «f + 0.5 we can see a change 
in the shape of the contours in N{CO). This change corresponds 
to an increase up to a factor of two in N{CO), with no notica- 
ble change in regions where logjo n > logjg «e + 0.5. From our 
fits of surface heating these regions have Fsurface < rmech/20. It 



^ The critical density «„- should not be confused with n^, defined ear- 
lier, which is the gas density where FnKch = Ftotai.surface 



is also noticed that the increase of N(CO) is more pronounced 
for high values of Go/n as Fmech increases (compare left and 
right panel of Fig.|6]l. We can explain this by looking at the the 
C^/C/CO transition zone, which gets shifted towards the surface 
of the cloud as Fmech dominates. This shift is more pronounced 
for higher Go, thus resulting in more CO in the column of H we 
are considering (A^(H) = 1.8 x lO^^cm"^). 

For conditions similar to that of the Milky Way, SFR ~ 
1 Mq yr"' which corresponds to a Fmech of 2 x 10"^" erg cm"^ 
s"' , N(CO) estimates would be marginally ( ~ 20% ) lower 
when Fmech is not taken into account. However, for starbursts 
with « ~ 10^ cm"^ and Go ~ 10^ and a significantly higher 
SFR (~ 50 Mq yr"') the A^(CO) estimates would be a factor of 
two less when the Fmech of 1 x 10"'^ erg cm"^ s"' is ignored 
(compare Fig.|7]to the panel corresponding to Fmech = 10"'** erg 
cm"^ s"' in Fig. lA.2l where the N{CO) contour is shifted to the 
right, from n ~ lO' to « ~ lO'*'^ cm""* ). For low metallicity 
grids we considered Z = 0.5Zq, which may represent systems 
such as dwarf galaxies (where the metallicity may be as low as 
0. 1 Z0),the increase in N{CO) may be up to one order of magni- 
tude, reaching 10'^ cm"- . This occurs even for a low mechanical 
heating rates corresponding to an SFR of ~ 1 Mq yr"'. The in- 
crease in N{CO) is visible in the central part of the parameter 
space (compare the grids when F^ech = 10"^' erg cm""* s"' and 
Tmech = 10"^" erg cm"' s"' in Fig. lA.4b . In contrast, for the high- 
est metallicity we considered, the Z = 2Z0 grids, which repre- 
sents a more likely condition for galactic centers, the increase is 
less prominent even for the extreme mechanical heating rates of 
Tmech = 10"" erg cm"-' s"' corresponding to ~ 500 Mq yr"' 
(see Fig. lA.3b . 

Thus, we conclude that we would potentially severely under- 
estimate the amount of CO in the system if mechanical heating 
is ignored, especially for metallicities less than Zq, i.e. in dwarf 
galaxies. 



3.2.2. H2O 

Until the launch of the Herschel Space Observatory, it was 
not possible to observe water lines in the far-infrared from 
extra-galactic sources, as their emission is blocked by the 
Earth's atmosphere. Now, water line emission has been de- 
tected by several Herschel key programs, and a striking exam- 
ple is pr ovided by the water li nes observed in Markarian 231 
(Gonzalez -Alfonso et al.L l2010l) . In this object, water lines in- 
tensities are comparable to those of the CO lines. This is very 
diff'er ent from what is o bserved for, e.g., the star-burst galaxy 
M82 dWeiB et al.L l2010h and the Orion Bar in the Milky Way 
(lHabartetal.Ll2010h" 

Water can be formed through a neutral-neutral reaction 
chain, O H- H2 -> OH + H, followed by OH -H H2 ^ H2O + H. 
In order to obtain significant abundances, the medium has to be 
molecular and sufficiently warm, as the reactions have barriers 
of ~ 2000 K. In X-ray exposed environments, it is possible to 
obtain a significant water abundance xh,o ~ 10"'° at high tem- 
peratures {T ~ 3000 K), produced by an ion-molecule reaction 
chain. A recent paper by Meijerink et al. (201 1) shows that water 
production can be enhanced by an order of magnitude in these 
warm atomic environments, when the grains are stripped from 
their ice layers, as is possibly the case in violent environments 
around an AGN. 

As we have seen in the previous sections, the temperatures 
in the molecular regions of the mechanically heated PDR can be 
increased to temperatures T » 100 K, which is favorable for 
the production of water through the neutral reaction chain. 
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Fig. 7. CO total column density up to A^(H) = 1 .8 x 10^^ cm - or 
Ay ~ 10 mag for Z - Zq without mechanical heating. 



In Fig. [8] the abundance of H2O is shown as a function of 
depth. The low-density models are marginally affected as long 
as Fmech < 10"^° erg cm"-' s"' . For higher mechanical heat- 
ing rates the H/H2 transition does not occur, thus H2O does not 
form. However, the impact on the abundance of H2O is more 
significant for the the high-density (M3 and M4), where the 



abundance of H2O is enhanced by at least one order of mag- 
nitude when Fmech > 10 '^ erg cm"^ s"' . The reference grid 
for the column density of H2O at solar metallicity (without me- 
chanical) heating is shown in Fig.|9] In comparing this refer- 
ence grid to Figs. IA.5l IA.6I and IA.7I (corresponding to grids for 
lower and higher metallicities), the column densities range be- 
tween A^(H20) ~ 10'"* (very low densities, n < 10' cm""* ) 
and up to 10'*'^ cm"-^ (n > 10^ cm"-' ). These column den- 
sities increase when the gas density if close to n — 11^ espe- 
cially for F^ech > 10"'^ erg cm"-' s"' where an optimum of 
(A^(H20) ~ 10'^ cm"^) is reached (see last row of Fig. lA.5) . 
When the mechanical heating rate becomes too high (in the re- 
gions where n « iic), the column densities drop again, because 
there is insufficient H2, which is destroyed through collisional 
dissociation, to form water Therefore, one finds that this opti- 
mum moves from lower to higher densities for increasing me- 
chanical heating rates. 

Thus, the high column densities of H2O accompanied by 
high temperatures make this molecule an excellent diagnostic 
tool for galaxy environments where mechanical heating domi- 
nates. 
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Fig. 6. CO abundance as a function of Ay for selected mod- 
els. In all the diagrams, the C/CO transition zone is shifted to- 
wards the surface as a function of increasing Fmech- For Ml 
and M2, the abundance of CO is below the displayed range for 
Fmech > 10"^*^ erg cm"-' s"' (these curves are not shown). The £. 
other curves almost overlap except for Fmech - 10"^*^ erg cm"^ ^ 
s~'. In models M3 and M4 the abundance of CO drops slightly 5 
beyond the transition zone compared to the case of no mechan- 
ical heating. The curves overlap except for Fmech = 10"'* erg 
cm"^^ s"' where the C^/C/CO transition does not fully take place. 
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3.2.3. HCN and HNC 

The critical density for the HCN and HNC 7=1-0 collisional 
excitation by H is ~ 1 x 10* and ~ 0.5 x 10* cm"-' respectively. 
Together with their relatively high abundance, makes them valu- 
able high gas density tracers. Here we will be first looking at the 
abundances of HNC and HCN, xhnc and xhcn respectively up 
to Ay = 20 mag. Then we will study the effect of mechanical 
heating on the column density ratios of these two species. 

Although included for completeness in Figs.[TO]and[TT] the 
low-density models Ml and M2 are of no physical interest, 
since they refer to densities at least two order of magnitude 
below the critical density of HCN and HNC. The more rele- 
vant high-density models M3 and M4 show (in contrast to the 
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Fig. 9. H2O total column density up to A^(H) - 1 .8x10" cm - or 
Ay ~ 10 mag for Z -Z^ without mechanical heating. 



low-density models), an increase by four orders of magnitude 
in jchcn when Fmech exceeds 10"'^ erg cm~^ s"'; whereas an 
increase of three orders in ahnc is obser ved. The incr e ase in 
the abundance of HCN w as also noted by iLoenen et al.l (l2008h 
and lMeiierink et all (1201 ll) . Mechanical heating causes the equi- 
librium temperature, in the molecular region, to increase from 
100 K, when r^ech = lO^'^erg cm^^ s"', to ~ 900 K when 
Tmech = 10"'* erg cm""* s"' (see Fig.|5]). Consequently HNC is 
converted more effectively to HCN via the reaction HNC + H 
— » HCN -H H whose activation barrier of ~ 200 K is surpassed 
dSchilke et al.L |1992|) . This explains the drop in the fractional 
abundance of HNC with respect to HCN with xhnc/^hcn ~ 0.1. 

The abundance of HCN and HNC increases by few orders of 
magnitude with increasing Fmech ■ This increase is because of the 
enhanced production of both species via the neutral-neutral reac- 
tion H2 + CN ^ HCN + H. Also HCN and HNC are formed with 
equal probability through reco mbinati o n reac tions of HCNH^, 
which has been also noted by lAaltol (l2005h . These two pro- 
cesses are more dominant than the reaction HNC H- H — > HCN 
H- H, which is just a re-shuffling the bonds of HNC that leaves 
■^HCN + -^^HNC constant. However HNC + H ^ HCN + H is re- 
sponsible for the bias in the abundance of HCN with respect to 
HNC. 

Now we look at the column density ratio A^(HNC)/A^(HCN) 
in order to gauge the effects of mechanical heating in the 
Go - n parameter space. These are computed by integrating 
xhnc and xhcn up to A^(H) = 1.8 x 10^^ cm"^ and taking 
their ratio. The same column density is considered for the other 
metallicities when the ratio A^(HNC)/A^(HCN) is computed. In 
Fig.[T2] the column density ratio of HNC to HCN is shown 
for the whole grid with no mechanical heating (pure PDR). 
Since HCN is a good tracer for den se molecular clouds with 



lO'* cm ^ (ISolomon et al.L Il992h . we focus on that part of 



n > 

the parameter space for this diagnostic. In the absence of me- 
chanical heating, in this range of density, A^(HNC)/A^(HCN) < 
1 for all Go. When mechanical heating is added, specifically for 
Fmech larger than 10"'^ erg cm""* s"', the column density ratio 
starts decreasing even more with almost no dependence on Gq. 
For instance, in the upper right comer A^(HNC)/A^(HCN) de- 
creases from 0.7 to 0.55 (see Fig. lA.8) . This decrease affects 
regions where log[o« - logio«c + 1-0. At logio« - logio«f 
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Fig. 10. HNC abundance as a function of Ay for selected models. 
In the top panels the curves for Fmech > 10~^^ erg cm~^ s"' are 
not shown since the C^/C/CO transition does not take place. 
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Fig. 11. HCN abundance as a function of Ay for selected models. 
The same convention of Fig.[TO]are used. 



(dashed vertical lines in Fig. lA.8l l F^ech = Fsurtace as described in 
Sect. l3.1.n so mechanical heating dominates everywhere in the 
cloud. Whereas at logjo « - logjg n^ + 1.0 mechanical heating is 
up to 100 times smaller than the total heating at the surface of 
the pure PDR case. For the largest F^ech considered 10"' erg 
cm"^ s"' the A^(HNC)/A^(HCN) decreases below 0.1 over the 



whole parameter space (see Fig. lA.8) . This is also observed for 
the other metallicities in Figs. lA.10l and lA.9l 
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Fig. 12. HNC/HCN integrated column density ratios (integrated 
up to Af(H) = 1.8 X \(f--cm-^ or Ay ~ 10 mag for Z = Zq) 
without mechanical heating. The column density ratio decreases 
from left to right. For dense molecular clouds, the ratio is ~0.7 
and it is almost constant as a function of Go. 



Thus, we find that the ratio A^(HNC)/A^(HCN) is sensitive 
to low mechanical heating rates, a factor of a hundred smaller 
than the total heating at the surface of a pure PDR. For higher 
mechanical heating rates that are 20 times smaller than those at 
the surface, A^(HNC)/A^(HCN) decrease at least by a factor of 
two. 



3.2.4. HCN/HCO+ 



'. This is caused by a slightly slower destruc- 
through the reaction HCO+ + e" ^ CO + 

10""^ erg cm""* 



ders of magnitude compared to the other reactions. However an 
increase of one to four orders of magnitude in xhco+ is observed 
for Ay < 5 mag. This contributes significantly to the column 
density of HCO^ in clouds with low A^(H), which could make a 
significant impact on the interpretation of line ratios for clouds 
with small column densities as opposed to those with high A^(H). 

Looking at the column density ratios of HCN to HCO^, 
A^(HCN)/A^(HCO+), we see in Fig. [14] a non-trivial dependence. 
The ratio varies from 10""^ (in the high-density and high-UV 
flux zone) to 10*^^ (in the high-density, moderate- to low-UV 
flux zone). 

An interesting feature of this column density ratio, quite 
different from the HNC/HCN situation, is the fact that the 
A^(HCN)/A?(HCO+) ratio increases to ~ 10 and ~ 10^ as 
opposed to a decrease of only a factor of two in the ratio 
A^(HNC)/A^(HCN) as a function of increasing Fmech ■ For high 
densities when F,j,ech > 10"^" erg cm"^ s~' the contours become 
almost parallel forlog[o«c < logjQ« < logjgnf + 1.0. This im- 
plies fliat the A^(HCN)/A^(HCO+) ratio does not depend on Go/n 
(as far as this ratio is concerned), but is dominated by mechan- 
ical heating down to values ~ 100 times smaller than the sur- 
face heating of a pure PDR. This is shown in Fig. lA. 11] where 
A^(HCN)/A^(HCO+) ranges between 10 and lO'' in the density 
range mentioned. Similar behavior is seen for the other metallic- 
ities considered. 

This, and the fact that A^(HCN)/A^(HCO+) is very sensitive 
to very small values of mechanical heating makes this ratio an 
excellent diagnostic for mechanically heated regions. 
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HCO^ is another high density gas tracer with yicr ~ 10^ cm"^ for — 
7=1-0 transition. It has been suggested that in AGNs, where 
X-ray heating dominates, HCN is less abundant than HCO^ 
Cepp & Dalgarno, 1996), suggesting that the HCN/HCO+ ra- 
tio is a diagnostic for the determination of the physical nature 
of a source object. A d etailed study of the effect of cosmic rays 
dMeiierink et al.Ll201 1]) shows that the abundance of HCO^ is in- 
hibited by high cosmic ray rates. Since cosmic rays have a large 
penetrating depth, they may raise the temperature of clouds up to 
large column densities. Since Fmech is added uniformly through- 
out the slab and the temperatures are likewise enhanced, a be- 
havior similar to that caused by cosmic rays should be seen. In 
a similar fashion as for HNC and HCN, we look at abundances 
and column density ratios with and without mechanical heating 
to see how it affects the A^(HCN)/A^(HCO+) ratios. 

In the low- density models (top row of Fig.[T3]l a slight in- 
crease in the abundance of HCO^ is observed up to Fmech - 
10"^^ erg cm"-' s"' ™- • 
tion rate of HCO^ 

H. For r,nech > 10""^ erg cm""* s"', the two orders of magni- 
tude drop in the abundance of HCO^ in the molecular region is 
caused mainly by ion-neutral reactions with H2O, HCN (which 
we demonstrated to increase in the previous section) and C. The 
abundance of HCO^ is unaffected by r^ech in M3 and M4 unless 
it is larger than 10"'^ erg cm""* s"' . For that amount of mechani- 
cal heating xhco+ becomes ~ 10"" beyond the radical region. In 
this case, the reactions H2O + HCO+ ^ CO + H3O+ and HCN + 
HCO"^ -^ HCNH^ -I- CO dominate (again because of enhanced 
abundances of H2O and HCN) in destroying HCO^ by two or- 
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4. Conclusion and discussion 

The goal of this paper is to assess when we should take into 
account mechanical heating in modeling the observed line emis- 
sion properties of the interstellar medium in PDRs. We have cal- 
culated the equilibrium thermal and chemical properties of in- 



11 



Kazandjian et. al: Diagnostics of PDRs with Mechanical Heating 




Fig. 14. HCN/HCO^ integrated column density ratios (inte- 
grated up to A?(H) = 1.8x10" cm-' or Ay ~ 10 mag for Z = Zq) 
without mechanical heating. 



terstellar clouds as ID semi-infinite slabs for a large range of 
hydrogen gas densities (1 < « < 10^ cm^S) and radiation fields 
(0.3 < Go < 10^). As a result, we covered more or less all con- 
ditions physically relevant for the ISM, from the diffuse ISM in 
the outskirts of quiescent galaxies (Go < 1, n = 1 - 100 cm"^), to 
the dense, UV-exposed environments (Go ~ 10^, « ~ 10^^ cm""*) 
of active galaxy centers. We have also explored three levels of 
heavy-element abundance, Z = 0.5, 1.0 and 2.OZ0, thus con- 
sidering the low metallicities of dwarf galaxies as well as the 
enhanced metallicities in galaxy centers. We have done this by 
adding a mechanical heating term resulting from dissipating tur- 
bulence of supemovae. The heating rates included are between 
10"""* < r„,ech < 10"'^ erg cm"-' s"', which corresponds to star- 
formation rates ranging from SFR ~ to ~ 1000 M© yr"' under 
the assumptions explained in section 2.1. We first explored the 
effect on the thermal balance at the surface and deep into the 
cloud at high column densities, and then the chemistry of a few 
commonly observed species in the ISM of galaxies, namely, CO, 
H2O, HCN, HNC, and HCO+. We demonstrated that mechani- 
cal heating radically changes the resulting chemical and thermal 
structure. The following conclusions can be drawn : 
(i) The equilibrium temperature in a PDR is significantly altered, 
i.e., increased by an order of magnitude, for ratios Fmech/n''^ ~ 
10"^^ erg cm""*-^ s"'. This is a lower bound in the sense that 
when Fmech/n''^ = 10""^ the temperature inside the cloud (even 
for the densest clouds n ~ 10^ cm""*) increases by a factor of 
10. Thus, Tmech mustht taken into account while solving for the 
equilibrium state of a PDR. In translating this relation into one 
relating gas density to the star formation rate, it can be expressed 
as SFR/n^^~ = 7 x 10"^. For higher values the temperature in- 
creases even more. Also, the H/H2 transition may not occur for 
very high mechanical heating rates (F,nech/n'^ 3> 10"^^), espe- 
cially for the low density clouds with n ~ 10^ cm"^. 
(ii) Mechanical heating rates that are two orders of magnitudes 
smaller than the surface heating nevertheless are high enough 
to increase the temperature by an order of magnitude in low- 
density clouds and up to a factor of two in high-density clouds 
n ~ lO^'cm"-'. For instance, inclusion of a mechanical heating 
rate of 10"'^ erg cm~^ s^' in the high-density cloud models M3 
and M4 of Meijerink & Spaans (2005) increases the tempera- 



tures inside the clouds from T ~ 30 K to ~ 100 K and from 
T ~ 70 K to ~ 110 K, respectively. Such heating rates Fmech 
are expected to occur in star-burst galaxies, and correspond to 
an SN rat e of approximately .3 per year and an SFR rate of 9.4 
Mq yr' (lLoenenetal.l(l2008h ). 

(iii) We observe that N{CO) is easily doubled when mechanical 
heating is included in modeling PDRs, especially for gas with 
« ~ lO-'cm"-'. As such, mechanical heating should be consid- 
ered, when interpreting observed CO data. We also observe that 
the implied column density strongly depends on metallicity. This 
column density decreases from 10'** cm"- at solar metallicity 
to 10'^ cm"^ at half-solar metallicity; whereas it increases up 



to 10 cm" at double-solar metallicity. This is particularly in- 
teresting for galaxy centers where the metallicity is expected to 
be higher than that of the solar neighborhood, 
(iv) Our results show that the brightness of emission lines ob- 
served in galaxy centers may at least partly be due to high 
mechanical heating rates. The models with Z = 2Z0 require 
higher column densities, and temperatures are enhanced by a 
factor 2 to 10. It also suggests that CO column densities in low- 
metallicity dwarf galaxies are even lower than assumed other- 
wise. Mechanical heating excites the gas to high temperatures, 
but the column density N{CO) is significantly lower (by almost 
an order of magnitude) compared to Z = Z©). In addition to that, 
the relative gas fraction is lower in dwarf galaxies, as supemovae 
easily disperse gas from the galaxy. 

(v) Like CO, H2O also exhibits an increase in column density. 
The temperature increase induced by mechanical heating results 
in implied column densities A^(H20) three orders of magnitude 
higher (up to ~ 10'^ cm"^) compared to models without mechan- 
ical heating. This gives H2O approximately the same column 
density as CO, for the highest mechanical heating rate. Based 
on this, H2O lines could also provide an excellent diagnostic 
in combination with CO, when studying mechanical heating in 
galaxies. 

(vi) HCN, HNC, HCO-H line ratios h ave been used t o dis- 
criminate be t ween PDRs and XDRs (iMeiierink et al.L 120071: 
iLoenen et all 120081) . Although we do not compute line intensity 
ratios in this paper, we have studied the signature of mechanical 
heating on the integrated column density ratios of HNC/HCN 
and HCN/HCO-H. We find that model HNC/HCN column den- 
sity ratios decrease by a factor of at least 5 for gas densities 
n > 10^ cm"^, when Fmech > 10"'^ erg cm""* s"'. This ef- 
fect is mitigated in clouds subjected to lower mechanical rates, 
but only at densities n < 10"* cm"^. HCN and HNC do not 
produce significant line emission at these densities as they are 
well below the critical density for excitation. On the other hand, 
the HCN/HCO^ column density ratio strongly depends on me- 
chanical heating. This ratio increases from ~ 10 to ~ 100, 
for dense molecular clouds (« > 10^ cm"^), when Fmech = 
10"'^ erg cm"-' s"' and up to lO'* for F^ech - 10"'^ erg cm"-' s"' 
for Z = 0.5, 1 .0, 2Zo. In general for these conditions the column 
densities of HCO^ decrease from ~ 10'"* cm"^ to ~ 10 
those for HCN increase from ~ lO'-^cm"- to 
whereas A^(HNC) drop to ~ 10'^ cm"- from ~ 10"*cm"^. The 
abundance of HCO^ drops by at least three orders of magnitude 
at column densities higher A^h > 1-8 x 10"' cm"- for the afore- 
mentioned Fmech- In galaxy centers, higher Av are easily possi- 
ble compared to Z = Z©. This translates into higher A^(H), which 
in turn would imply a smaller A^(HCO^) compared to A^(HCN). 
Hence resulting in even higher HCN/HCO+ column density ra- 
tio. The lack of HCO^ in combination with bright HCN line 
emission might be useful in tracing regions with high mechani- 
cal heating rates. 
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The integrated column density ratios do not translate di- 
rectly into line intensity ratios. Instead, one should compute 
these ratios by solving the radiative transfer equations. We will 
consider predictions for the actual line intensities of HCN, 
HNC and HCO^ in the next paper, where we will also as- 
sess time-dependent effects. The equilibrium treatment is cer- 
tainly valid for the low-metallicity dwarf galaxies, and the 
benchmark-metallicity of solar neighborhood galaxy discs, but 
this may not be the case for the high-metallicity galaxy cen- 
ters, where the chemical, thermal, and dynamical timescales are 
quite comparable in the inner < 100 pc, especially at densi- 
ties w < 10^ cm~^. For examp le, the formation timescale of H2 
(iHollenbach & TielensLll999l) is < 5 Myr at an ambient gas den- 
sity n ~ 10^ cm -^ and a kinetic temperature ofT~ 100 K. The 
stellar orbital periods (which stir up the ISM) are comparable to 
or shorter than this, namely between 0.1 and 15 Myr Hence, a 
time-dependent treatment is necessary, and we aim to accom- 
plish this in upcoming papers by adapting the current PDR- 
XDR code. We will allow for time-dependent chemical evolu- 
tion, and couple it with an SPH code through AMUSE. As such, 
the chemistry, gas and stellar dynamics are being evolved self- 
consistently. 

Acknowledgements. M.V.K would like to thank A.G.G.M. Tielens, M. 
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References 

Aalto, S. 2005, in lAU Symposium, Vol. 231, Astrochemistry: Recent Successes 

and Current Challenges, ed. D. C. Lis, G. A. Blake, & E. Herbst, 261-270 
Bakes, E. L. O. & Tielens, A. G. G. M. 1994, ApJ, 427, 822 
Bayet, E., Williams, D. A., Hartquist, T. W., & Viti, S. 201 1, MNRAS, 414, 1583 
Black, J. H. & Dalgamo, A. 1977, ApJS, 34, 405 
Bohlin, R. C, Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132 
Dahlen, T. & Fransson, C. 1999, A&A, 350, 349 
de Jong, T., Boland, W., & Dalgamo, A. 1980, A&A, 91, 68 
Gonzalez-Alfonso, E., Fischer, J., Isaak, K., et al. 2010, A&A, 518, L43 
Habart. E., Dartois, E., Abergel, A., et al. 2010, A&A, 518, L116 
Habing, H. J. 1969, Bull. Astron. Inst. Netherlands, 20, 177 
HoUenbach, D. J. & Tielens, A. G. G. M. 1999, Reviews of Modem Physics, 71, 

173 
Israel, F. P. 2009, A&A, 506, 689 
Kaufman, M. J., Wolfire, M. G., Hollenbach, D. J., & Luhman, M. L. 1999, ApJ, 

527, 795 
Kroupa. P 2001, MNRAS, 322, 231 

Leitherer, C, Schaerer, D., Goldader, J. D., et al. 1999, ApJS, 123, 3 
Lepp, S. & Dalgamo, A. 1996, A&A. 306. L21 

Loenen, A. F., Spaans, M., Baan, W. A., & Meijerink, R. 2008, A&A, 488, L5 
Maloney, P R., Hollenbach, D. J., & Tielens, A. G. G. M. 1996, ApJ, 466, 561 
Meijerink, R. & Spaans, M. 2005, A&A, 436, 397 
Meijerink, R., Spaans, M., & Israel, F P 2007, A&A, 461, 793 
Meijerink, R., Spaans, M., Loenen, A. F, & van der Werf, P P 201 1, A&A, 525, 

A119+ 
Papadopoulos, P P, van der Werf, P, Isaak, K., & Xilouris, E. M. 2010, ApJ, 

715,775 
Pelupessy, F. I. 2005, PhD thesis, Leiden Observatory, Leiden University, 

PO. Box 9513. 2300 RA Leiden, The Netherlands 
Pelupessy, F I. & Portegies Zwart, S. 2012, MNRAS, 420, 1503 
Perez-Beaupuits, J. P, Wada, K., & Spaans, M. 201 1, ApJ, 730, 48 
Portegies Zwart. S., McMillan, S., Harfst, S., et al. 2009, New A, 14. 369 
Rollig, M. 201 1, A&A, 530, A9+ 
RoUig, M., Ossenkopf, V, Jeyakumar, S., Stutzki, J., & Sternberg, A. 2006, 

A&A. 451,917 
Schilke, P, Walmsley, C. M., Pineau Des Forets, G., et al. 1992, A&A, 256, 595 
Solomon, P M., Downes, D., & Radford, S. J. E. 1992, ApJ, 387, L55 
Tielens, A. G. G. M. & Hagen, W. 1982, A&A, 1 14, 245 
Tielens, A. G. G. M. & Hollenbach, D. 1985, ApJ, 291, 722 
Wada, K. & Nomian, C. A. 2002, ApJ, 566, L21 
Wada, K. & Tomisaka, K. 2005, ApJ, 619, 93 
WeiB, A., Requena-Torres, M. A., Gusten, R., et al. 2010, A&A, 521, LI 



13 



Kazandjian et. al: Diagnostics of PDRs with Mechanical Heating 



r„=io-^' 



r„=io' 



r„=io-^' 



5 
4 

a 

" 3 

D> 
O 

2 
1 


5 

4 

o° 

" 3 

01 

o 

2 

1 



2 



1 




r„=io 




r„=io 



r..=io' 



r..=io' 



r„=io' 



r..=io' 



"T 1 rn 



2 3 4 5 
log^Jn/cm ] 



2 3 4 
log,„[n/cm ] 





1 
1 

1 


.1 


Si 

1 


1 


1 1 1 
1 1 1 


1 

S 

g 


s 


7 


III, 
III. 


J 

1 


/ 

a 



2 3 4 5 
log,„[n/cm ] 



^-24 



Fig. A.l. Temperature contours at Ay - 20 for different values of mechanical heating. As mechanical heating is increased from 10" 
to 10"'*, the temperatures increase from left to right in the grid as function of increasing mechanical heating. The temperatures to 
the left of the highest temperature contour line are higher than the temperature of the contour curve itself, which are usually larger 
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Fig. A.2. logjQ CO total column density (integrated up to A^(H) = 1.8 x 10 cm or Ay ~ lOmag for Z — Z©) for different values of 
mechanical heating. 
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Fig.A.4. logio CO ^o^^ column density (integrated up to A^(H) = 1.8 x 10 cm or Ay ~ 5mag for Z — O.SZg) as a function of 
mechanical heating for the low metallicity case. 
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Fig. A.6. logjQ H2O total column density (integrated up to A^(H) = 1.8 x lO^^cm ^ or Ay ~ 20mag for Z = 2Zq) as a function of 
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Fig. A.7. logio H2O total column density (integrated up to A^(H) = 1 .8 x 10 cm or Ay ~ 5mag for Z = O.SZq) as a function of 
mechanical heating for the low metallicity case. 
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Fig. A.8. HNC/HCN integrated column density ratios (integrated up to A^(H) = 1.8 x lO^^cm ^ or Ay ~ lOmag for Z = Zq) as a 
function of mechanical heating. 
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Fig.A.9. HNC/HCN column density ratios (integrated up to A^(H) = 1.8 x 10~^cm ~ ox Ay ~ 20mag forZ - 2Zq) as a function of 
mechanical heating for the double solar metallicity case. 
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Fig. A.IO. HNC/HCN integrated column density ratios (integrated up to A^(H) = 1.8 x lO'^cm ^ or Ay ~ 5mag for Z - O.SZq) as a 
function of mechanical heating for the low metallicity case. 
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Fig. A.ll. logjo HCN/HCO^ integrated column density ratios (integrated up to A^(H) = 1 .8 x 10 cm ^ or Ay ~ lOmag for Z = Zq) 
as a function of mechanical heating for solar metallicity. 
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Fig. A.12. logjo HCN/HCO^ column density ratios (integrated up to A^(H) 
function of mechanical heating for the double solar metallicity case. 
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Fig.A.13.1og[oHCN/HCO+ integrated column density ratios (integrated up to A^(H) - l.SxlO^^cm -orAy ~ SmagforZ = O.SZq) 
as a function of mechanical heating for the low metallicity case. 
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